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Abstract 

We report measurements of the branching fraction and CP violation parameters in B° — ► D + D~ 
decays. The results are based on a data sample that contains 535 x 10 6 BB pairs collected at the 
T(4S) resonance, with the Belle detector at the KEKB asymmetric-energy e + e~ collider. We obtain 
[1.97 ± 0.20 (stat) ± 0.20 (syst) ] x 10" 4 for the branching fraction of B° -► D+D~. The measured 
values of the CP violation parameters are: S = -1.13 ± 0.37 ± 0.09, A = 0.91 ± 0.23 ± 0.06, 
where the first error is statistical and the second is systematic. We find evidence of CP violation 
in B° — > D + D~ at the 4.1 a confidence level. While the value of S is consistent with expectations 
from other measurements, the value of the parameter A favors large direct CP violation at the 
3.2c7 confidence level, in contradiction to Standard Model expectations. 

PACS numbers: 13.25.Hw, 11.30.Er 
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Within the Standard Model (SM), CP violation (CPV) arises from a complex phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix V [1]. The dominant contri- 
bution to B° — > D + D~ decays is the tree-level b — > ccd transition shown in Fig. QJa). If this 
diagram is the only contribution, then the mixing- induced CPV parameter for B° — ► D + D~ 
is — sin 20x, where 4>i — ar g[ — VcdV* b /V t dV t l\, while the direct CPV term A is zero. The 
second-order gluonic penguin contribution, shown in Fig. QJb), is expected to change the 
value of the parameter S by less than a few percent and increase the value of A to about 
3% [3]. However, particles from physics beyond the SM may give additional contributions 





(a) 



(b) 



FIG. 1: The tree (a) and the penguin (b) contributions to the B° — > D + D decay. 



within the loop diagrams mediating flavor- changing b — > d transitions. Such contributions 
may potentially induce large deviations from the SM expectation for time-dependent CP 
asymmetries. As sin 2<fii has already been determined with high precision by measurements 
in b — > ccs charmonium modes j^, the objective here is to focus on deviations from ex- 
pectations in b — > ccd transitions. Similar studies have been carried out for B° — > D^D"^ 
decays, which involve the same quark level weak decay @, Q, 0, l^]. 

The CPV parameters S and A can be measured from the At distribution of B° — > D + D~ 
decays, 



V 

' s: 



-|At|/r 



At 



-{1 + g[5sin(AmAt) + ^lcos(AmAt)]} 



(1) 



where At = tcp — ttag is the time difference between decays of the two B mesons arising from 
the T(4f>). The parameters tcp and t tag are the proper decay times of the corresponding 
B mesons, r is the B° meson lifetime and Am is the mass difference of the two B mass 
estates Q. The flavor , is denned flea the flaa. state of the ta gg ing B meson: 
g = +1 (-1) for B tag = 5° (5°). 

The results presented here are based on a data sample that contains (535 ± 7) x 10 6 
BB pairs, collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 
8 GeV) collider KEKB operates at the T(45) resonance (y/s = 10.58 GeV) with a 
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peak luminosity that exceeds 1.7 x 10 34 cm~ 2 s -1 . At KEKB, the T(4S) is produced with a 
Lorentz boost of f3j = 0.425 nearly along the electron beam line {—z direction). Since the 
B° and B° mesons are approximately at rest in the T(4S) center-of-mass (CM) system, At 
can be determined from the displacement in z between the Bqp and B taig decay vertices: 
At ~ (z C p - zt ag )/(3yc = Az/Pyc. 

The Belle detector [l2J is a large-solid-angle magnetic spectrometer that consists of a 
silicon vertex detector, a 50- layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons. Two 
inner detector configurations were used; a 2.0 cm radius beam pipe and a 3-layer silicon 
vertex detector was used for the first 152 x 10 6 BB pairs and a 1.5 cm beam pipe, a 4- 
layer silicon detector and a small-cell inner drift chamber were employed for the remaining 
383 x 10 6 BB pairs (l3|. 

D mesons are reconstructed using the D + — > K~7t + tt + and D + — > Kstt + decay modes 



14| . In this paper, the shorter notation Kmc is used when both D mesons are reconstructed 



in the Kirn channel while Ks7r is used when at least one of the D mesons is reconstructed 
in the K„ ehanne,. Charged tracks that ate not pos.t.ve.y identified as electrons Q 
and satisfy a loose requirement on the impact parameter relative to the interaction point 
(IP) are considered as pion and kaon candidates. For charged particle identification (PID), 
we combine information from the CDC, TOF and ACC counters into a likelihood ratio 
C(K ) /^(K^ + £(7r ± )]_ I _which is required to be greater than 0.55 for kaon and less than 



0.9 for pion candidates 16|. K$ candidates are reconstructed in the K$ — > 7r + 7r decay 
mode; the pion combination is required to have an invariant mass within 30MeV/c 2 of the 
nominal Kg mass and a vertex displaced from the IP. The mass of the meson candidate 
is required to be within 10MeV/c 2 (2.4er) of the nominal D ± mass. We select B meson 
candidates using the energy difference AE = E* B — E^ eam and the beam-energy-constrained 
mass M bc = \J (E^ cam /c 2 ) 2 - (p* B /cf, where E* B , E^ eSLm , and p* B are the B meson energy, 
the beam energy, and the B meson momentum, respectively, in the CM system. 

The Ks decay vertex is fitted from two pion tracks. The D + meson decay vertex is 
fitted from three charged tracks or from the Ks and 7r + track. The mass of the K^-k + -k + 
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or Kstt + combination is constrained to the D + meson mass to obtain better M^ c and AE 
resolutions. The B° decay vertex is reconstructed from the two D meson tracks and the IP 
information. All remaining charged tracks are used to determine the decay vertex of the 
tag-side B meson. A loose requirement on the quality of the vertex fit is applied for both 
B mesons. The reconstruction of the B tSug vertex, vertex quality and flavor tagging are not 
required for the branching fraction measurement. 

The flavor of the accompanying B meson is determined from its decay products. Events 
are divided into six r-bins according to the tagging quality r. The value of r ranges from 
for events with no flavor information to 1 for unambiguous flavor assignment. Due to the 
imperfect flavor tagging, the distribution P sig of Eq. flTJ is modified to 

e -|At|/r 

V sig = — {1 - qAw + q(l - 2w) (2) 

4r 

[S sin( Am At) + A cos( Am At)] } , 

where w is the wrong tag fraction, and Aw is the difference between the wrong tag fractions 
if the -Stag meson is a B° or B°. The values of w and Aw for each of the six bins in the 
sagging quality parameter r are determined separately using flavor specific B meson decays 
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Continuum events are suppressed by forming a likelihood ratio from cosds, where 6b 
is the polar angle between the B meson direction in the CM system and the beam axis, 
and a variable based on a combination of sixteen modified Fox- Wolfram moments with the 



scalar sum of transverse momentum [18j. Note that since the BB and continuum events 
have significantly different distributions in the tagging quality variable r, the continuum 
suppression cut varies for events in different r-bins. 

After applying all of the event selection criteria, 16% of the signal events have more than 
one B° candidate. The B° with the smallest value of (Am D + /a D +) 2 + (Am D - /a D -) 2 is 
selected as the best candidate, where Amp = M Knn / Ks7T — m D is the difference from the 
nominal D meson mass and cr D ± are the widths of the signal peak in the Mx n n/K S TT mass 
distribution. 

The signal yield is obtained from an extended unbinned 2D maximum likelihood (ML) 
fit of the M bc and AE distributions in the range M bc > 5.20GeV/c 2 and -0.05 GeV < 



AE < 0.10 GeV. A Gaussian function for the signal and an ARGUS 19] function for the 



background are used to describe the Mb c distribution. For the parameterization of the AE 



distribution we used two Gaussians with the same mean value to describe the signal and 
a linear function to describe the background. The fraction and the width of the wider 
Gaussian were fixed to the values obtained from Monte Carlo (MC) simulated signal decays 
Q. Non-resonant B° -> D-K°n + and B° -> D'T{* (S92)tt + decays are found to be a 
possible source of background peaking in the M^ c and AE distributions. The amount of 
this background was estimated from the D + mass sidebands in data and subtracted from 
the signal. We estimate the number of non-resonant decays in the signal region (N nr ) to be 
2.0 ± 1.8 and 1.4 ± 1.0 for the Kirn and Kgit channels, respectively. The fit yields 150 ± 15 
events in the signal peak, where the error is statistical only. The Mb c and AE distributions 
of reconstructed events and the projection of the fit result are shown in Fig. [2J The signal 
yields from separate fits to the Ktttt and Ksft decay modes are given in Table [B 
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FIG. 2: Distributions for the reconstructed events in M^ c c 2 (a) and AE (b). The full (dashed) 
curves show the projections of the result of the 2D unbinned maximum likelihood fit for all (back- 
ground) events. 



The combined branching fraction is calculated from the total number of reconstructed 
events and the average reconstruction efficiency, and is found to be B(B° — > D + D~) = 
1.97 ± 0.20 (stat) ± 0.20 (syst) ] x 10~ 4 , which is consistent with previous measurements 
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22j and has better accuracy. The uncertainty in the D meson branching fractions results 



in a 5% systematic error. The error in the pion and kaon track reconstruction efficiency was 
estimated using partially reconstructed D* decays. The errors are added linearly for all six 
pion and kaon tracks, which yields a 6% uncertainty. The difference in PID efficiency for 
the simulated and real data is approximately 1% per track, which gives a 6% uncertainty. 



TABLE I: The product of D branching fractions Bj^+ x Bp-, the detection efficiency e, the number 
of events in the signal peak iVp ea k and the expected amount of the combinatorial background Ab cg 
in the 5.27GeV/c 2 < M bc < 5.29GeV/c 2 and \AE\ < 0.03 GeV region, as extrapolated from the 

fit. 



channel 


B D + x B D - 




Apeak Abcg 


Kirn 


(0.904±0.065)% 


12.6 


124.1±13.6 110.8±2.6 




(0.204±0.015)% 


12.1 


25.7± 5.7 13.8±0.9 



Smaller contributions come from the uncertainty in the K$ selection efficiency (1%), the 
number of BB events (1.3%) and the number of non-resonant decays (1.5%). The total 
systematic error of 10% is obtained from the quadratic sum of these uncertainties. 

Time-dependent CP violation parameters are determined by an unbinned ML fit to the 
At distribution of 219 events in the signal region 5.27GeV/c 2 < M bc < 5.29GeV/c 2 and 
|Ai?| < 0.03 GeV. The At distribution for signal events V S i g described by Eq. ([2]) is modified 
by the inclusion of the background contribution and resolution effects. The event-by-event 
likelihood is given by 

£cv = /sig Kig ® Tl + /nr v m ® 71 + f hcg V hcg ® 7Z hcg . (3) 

Subscripts sig, nr and beg refer to signal, non-resonant and combinatorial background com- 
ponents, respectively. The fractions /; = /i(Mb c , AE, r) are determined on an event-by-event 
basis, / sig + / nr + /beg = 1. The function 71 describes the detector resolution of the At mea- 
surement. It takes into account the error in the determination of both B meson vertices as 
well as an additional kinematic smearing due to the momentum of the B meson in the CM 
system and the smearing of the tag-side vertex due to the tracks originating from the sec- 
ondary vertices. An additional wide Gaussian component with o w 20 ps is added to describe 
a small fraction of events (about 1 %) with poorly reconstructed vertices. A more detailed 



description of the resolution function parameterization can be found in Ref. [23|. Resolution 
parameters for the Bcp meson vertex are determined from a fit to the At distribution of 
kinematically similar B° D+D~ decays. 

The fraction of the non-resonant decays / nr is assumed to be proportional to the signal 
fraction, / nr = a/ sig , where a = N nr /(N pea _ k — N nr ) and a KlTn = 0.016, = 0.059. The 
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At distribution of the non-resonant B° — > D~K tt + or B° — > D~ K*° (892)tt + background is 
described by the At distribution for signal with the parameters S and A set to zero. We 
include the effect of possible CP asymmetry of these modes in the systematic error. About 
half of the combinatorial background events come from BB decays (b — > c transition), 
which have an exponential decay At distribution. The other half are continuum e + e~ — ► qq 
events, for which the At distribution contains a 5-function component. Therefore, the At 
distribution of the background is described by 



beg 



p-|At|/Tbcg 



2r b 



eg 



(4) 



The background resolution function TZ^cg is taken to be a Gaussian. Parameters describing 
the background distribution are obtained from a fit to the At distribution of the data 
sideband, M bc < 5.27GeV/c 2 and AE > 0.06 GeV. 

From an unbinned fit to the measured At distribution described by Eq. ([3]), we obtain 
the CP violation parameters for B° — * D + D~ , 



S = -1.13 ±0.37 ±0.09 and 
A = +0.91 ±0.23 ±0.06, 



(5) 



where the first error is statistical and the second is systematic. The At distributions are 
shown in Fig. [3j The main contributions to the systematic error are fit bias (0.06 for S and 




FIG. 3: The At distribution for events with good tagging information (r > 0.5) when the tag-side 
S-meson is reconstructed as B° (a) or B° (b). The full and dashed curves show the projection of 
the fit result and background contribution, respectively. 



0.02 for ^4), uncertainties in the resolution function (0.04 for S and 0.03 for ,4) and signal 
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fraction (0.035 for S and 0.015 for .4). Other uncertainties come from the errors on the 
parameters r and Am (0.023 for S and 0.007 for A), wrong tag fractions (0.017 for S and 
0.014 for ^4), description of background At distribution (0.01 for S and ^4), fraction and 
possible CP asymmetry of the non-resonant background (0.02 for S and 0.03 for .4), the 



effect of tag-side interference 



24j (0.01 for S and 0.03 for .4) and requirements on the vertex 



quality and the fitting range (less than 0.01 for S and 0.01 for ^4). 

To test the consistency of the fitting procedure, the same analysis was applied to the 
B° — > D+D~ control sample. Since there is only one decay amplitude at the tree level and 
the leading penguin contributions have the same CKM structure as the tree contribution, no 
CPV is expected for this decay. The result is consistent with no CPV, S = —0.064 ± 0.094 
and A = 0.091 ± 0.060, where the error is statistical only. We also fit the background 
sample (Mb c < 5.27GeV/c 2 and AE > 0.06 GeV) for a possible CP asymmetry and find 
none: A = —0.01 ± 0.06 and S = 0.03 ± 0.10. In addition, a time-integrated fit for the 
parameter A was performed to validate the result in B° — > D + D~ decays. The signal yield 
was determined separately for events tagged as B tSug = E>° and B tag = B° for each of the six 
r-bins. The fit yields A = 0.86 ± 0.32, which is consistent with the time-dependent result. 



We use the Feldman- Cousins frequentist approach [25j to determine the statistical sig- 
nificance of our measurement. In order to form confidence intervals, we use the A and S 
distributions of the results of fits to the MC pseudo-experiments for various input values of 
A and S in a similar way as described in Ref. 26j. The systematic errors and possibility 
of tails that are wider than Gaussian tails are taken into account. The case of no CPV, 
S = A = 0, is ruled out at the 4.1cx confidence level (CL). The case of no direct CPV, 
A — 0, is excluded at more than 3.2a CL for any value of the parameter S. 

In summary, we measure the branching fraction for B° — > D + D~ decays to be 
(1.97 ± 0.20 ± 0.20 ) x 10~ 4 , superseding our previous measurement [2l|. We obtain values 
for the CP parameters S = -1.13 ±0.37 ±0.09 and A = 0.91 ±0.23 ±0.06 and rule out the 
CP-conserving case, S — A — 0, at the 4.1 a confidence level. The value of S is consistent 
with measurements of b — >• ccs modes |ld |. In addition, we observe evidence for direct CP 
violation at the 3.2 a confidence level. Some extensions of the SM predict large contributions 
to the CP violating phases in b — > ccd decays that are consistent with our result 27]. Our 

n 

measurement differs from a previous measurement by the BaBar collaboration |8J by about 
2.2 a. 
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